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Abstract 

I consider the case of two interacting scalar fields, and if), and use 
the path integral formalism in order to treat the first classically and the 
second quantum-mechanically. I derive the Feynman rules and the resulting 
equation of motion for the classical field which should be an improvement of 
the usual semi-classical procedure. As an application I use this method in 
order to enforce Gauss's law as a classical equation in a non-abelian gauge 
theory. I argue that the theory is renormalizable and equivalent to the usual 
Yang-Mills as far as the gauge field terms are concerned. There are additional 
terms in the effective action that depend on the Lagrange multiplier field A 
that is used to enforce the constraint. These terms and their relation to the 
confining properties of the theory are discussed. 



1 Introduction 



Consider the case of two scalar fields, and ip, interacting with the action 



L = \{d<pf - - + \{d^f - - - (2) 



If one wishes to treat ip quantum-mechanically the usual, semi-classical, pro- 
cedure leads one to replace all terms containing ip in (3) by their quantum 
expectation values. There are, however, several problems associated with 
this procedure, conceptual as well as technical: One may use a background 
field method in order to derive the first (one-loop) quantum corrections to 
the effective potential and associated mass terms. If, however, a higher or- 
der calculation is needed, one would have to modify the propagator using 
the one-loop result, and be careful in order to avoid double-counting for any 
relevant diagrams. The procedure depends on the model examined and on 
relative order-of-magnitude estimates. One would like to have a more self- 
consistent approach in order to incorporate quantum effects on classical fields 
and vice-versa. 

Here I use the path integral formalism in order to describe this problem 
and I get the resulting Feynman rules that enable one to study the interac- 
tions between the classical and quantum fields. The result is what one would 
expect: the classical field propagates only in tree diagrams, not in loops, and 
the quantum field propagates as usual, providing the quantum corrections. 
This leads to an effective equation of motion for the classical field that should 
be an improvement upon the semi-classical procedure. The Feynman rules 
presented here reorganize the entire perturbation series and enable one to 
treat these problems in a self-consistent manner that can easily be extended 
in other models. 




(1) 



where 



The classical equation for <j) is obtained from 




(3) 
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The path integral approach has been used before in order to treat pure 
classical mechanics [1, 2] and investigate problems of classical behavior in 
quantum field theory [3, 4]. The formalism developed here has many similar- 
ities to these previous works. The main method is an extension of [1, 2] to the 
case of interactions between classical and quantum fields (there are, however, 
some important differences even in the case of purely classical fields). 

In Sec. 2 I develop the main formalism, derive the Feynman rules and 
get an effective action from which the equation of motion for the classical 
field can be obtained. An important advantage of this method is that it is a 
self-consistent procedure that can be used in order to calculate higher order 
effects. 

In Sec. 3 I use this method in order to treat Gauss's law as a classical 
equation in a non-abelian gauge theory. Because of the asymmetry of the 
Feynman rules described here, the effective action contains, besides the usual 
Yang-Mills terms, additional terms that depend on the Lagrange multiplier 
field that is used to enforce the classical constraint. These terms are of 
the Coleman- Weinberg type [8], and are reminiscent of phenomenological 
effective actions that have been used in order to describe confinement [11, 
12, 13]. They have similar interpretation here, where one can also see the 
region of validity of perturbation theory. 

Since the method presented here is new I include two Sections of com- 
ments. In Sec. 4 I discuss the applications to non-abelian gauge theory and 
in Sec. 5 I make some general comments regarding this method. 



2 Path integral and effective action 

In order to calculate 

Z(J, J') = J [d(f>] [#]<S(0 - cl ) exp (i J L(<j), V) + J0 + JV) , (4) 

where c i is the solution of (3), I use the Lagrange multiplier A and ghost 
fields c, c, similarly to the work in [1, 2], and adding another source A, the 
path integral to be evaluated becomes 

Z(J, J', A) = / [d<t>] [#] [dX] [dc] [dc] exp (i [ L + J<j> + J'i/j + AA) (5) 
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where 




(6) 



is the modified Lagrangian with the corresponding modified action S. For 
the simplest case of two interacting scalar fields with (2) we get 



where K — — (□ + m 2 ). The propagators and the vertices can be deduced 
from here. For the A and fields we get 



where N is a normalization factor, independent of the sources, and G = 
l/(k 2 — m 2 + ie), in momentum space, is the usual Feynman propagator. 
Accordingly there is no — propagator, there are, however, a mixed A — 
propagator equal to G, and a A — A propagator equal to —G. The remaining 
ip — ip and ghost propagators as well as the various vertices are as usual from 



One can now check: at one loop order the loops with the A — propagator 
cancel with the ghost loops and similar cancelations exist in higher loops, 
loops with the A — A propagator do not appear because of the Feynman rules 
of the modified action (the vertices are at most linear in A) with the final 
result that the field does not have quantum corrections but only propagates 
classically through tree diagrams. A typical line with the A and fields is 
either the sum of two A — and one A — A propagator, or a single A — 
propagator, in both cases equal to G. The ip field, of course, propagates also 
in loops as a genuine quantum field and gives the quantum corrections to the 
classical field 0. 

I should note here that the propagator G that we get for the classical field 
is the Feynman propagator, and not the retarded one that is usually employed 
in classical mechanics. There are two reasons for that: first, the boundary 
conditions used in the path integral (the field goes to zero at infinity) are 
different than the ones usually employed in classical mechanics (the field 
configuration is given at an initial time). One can check with a more careful 
evaluation of (8) that we get, indeed, the Feynman prescription. This is true 



S = S+ f (\K(I) + cKc - |rA0 3 - £'A0^ 2 + |c0 2 c + g'c^ 2 c) , (7) 



Jx \ 





(7). 
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even if we have only the classical field in our theory. One can presumably 
use the path integral formalism with different boundary conditions in order 
to attack purely classical problems. Then the retarded propagator would 
probably emerge, as in [2]. A second, physical reason that is relevant here, 
is that we want to study interactions between the classical and quantum 
fields. The possibility of particle creation and annihilation is essential for 
both the classical and quantum fields. Our classical field, therefore, admits 
both particles and antiparticles propagating classically. 

Renormalization of the theory proceeds as usual. All divergent terms 
come from loops of the quantum field if) and depend on the quantum-classical 
coupling g' or the quantum coupling g", not on g. The classical parameters g 
and m get also renormalized. However, pairs of terms like g\<f> 3 and g4> A in the 
modified action S have the same divergencies associated with, accordingly 
the renormalization procedure does not affect the classical nature of the field 



I will now proceed to show how the classical equation for gets modified 
in the presence of quantum interactions. If we use the generating functional 
Z(J, J', A) to construct W(J, J', A), the generating functional for connected 
diagrams, and from that r(0, ip, A), the effective action, we see that the 
appropriate equation is 



This can be verified if we consider the theory with only the classical field 
and use the properties of the Legendre transformation in order to express 
the classical equation SS/5(f) = —A = in terms of the effective action. So 
the equation of motion for the classical field <f> can be taken from the tadpole 
one-particle irreducible graphs with one external field A that contain if) loops 
and lines of the A and fields, but not loops of the classical field. For the 
case of the two interacting scalar fields of (2) we get at one loop order 




(9) 





(10) 



d<f) 



where 



t/ 9 ,a , 9' ,2 ,2 , * f d4k , 



( 



( fc 2_ m 2 )(A;2 _ M 2 ) + 

(k 2 - ml) 



) 



(11) 
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is the usual one loop effective potential of the original theory without the <p 
loop (mj = m 2 + \g<p 2 + g'tp 2 , M} = M 2 + g'<p 2 + \g"^ 2 ). In fact (11) also 
contains the loop with only the ip field. Since this is, however, ^-independent 
it does not contribute to (10). One can check that when the quantum-classical 
coupling g' is zero this reduces to the ordinary equation for a classical Klein- 
Gordon field. Regularization and renormalization are performed as usual, as 
was discussed above. 

This result can, of course, be also obtained using a background field 
perturbation theory. However, if one wishes to study higher order corrections, 
perturbation theory has to be reorganized at any order so as to incorporate 
the previous order result and avoid any multiple-counting. This problem does 
not appear here since we have a well-defined effective action, that depends 
on the auxiliary field A, with fixed Feynman rules from the beginning. 

The resulting equation for the classical field is written as an equation 
of motion from an effective action that contains what one would expect: 
arbitrary loops of the quantum field but not loops of the classical field. It 
depends on the vacuum expectation value of the quantum field ip, enabling 
one to also study problems of symmetry breaking. Even for ip — 0, of course, 
it provides the quantum corrections to the classical equations of motion. 
What is more important, this method allows one to include higher order 
corrections self-consistent ly through the effective action formalism. 

The complete effective action also describes the quantum properties of 
the field ip. One can then determine the effects of the classical field on the 
quantum field via its effective potential or other terms. There are also higher 
order terms that involve powers or derivatives of the auxiliary field A. Their 
relevance, if any, to the combined dynamics of the system is not clear from 
this work. In the simple example described here, since there is no symmetry 
breaking, we have to set A = anyway in order to derive the effective equation 
of motion. In cases with symmetry breaking, however, these terms may turn 
out to be important. One interesting case will be described in the next 
Section in the context of the non-abelian gauge theory, where they may be 
relevant to the confining properties of the theory. 
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3 An application in non-abelian gauge theory 



The formalism of the preceding section applies strictly in the case of two 
interacting fields. Here, however, I will show how it can be used in the 
case of the non-abelian gauge theory, in order to treat Gauss's law as a 
classical equation. The strategy will be the following: I will choose first a non- 
covariant gauge in which pure Yang-Mills is well-defined and renormalizable, 
and use the Feynman rules derived here in order to map a sector of the theory 
onto the usual Yang-Mills, and identify the remaining (or missing) terms as 
additional contributions to the effective action, that depend on A. 
For the non-abelian gauge theory with action S and Lagrangian 

L = --F°:F atxv (12) 

where F* v = d^Al — d u A a ^ + gf abc A b ^A c v , I will use Lagrange multipliers A" 
and ghost fields c a , c a , and get the modified action 

s - s+ h'^i§K cb (13) 



in order to treat Gauss's law 

5S 



5A 



= -D ■ E = (14) 



classically (S and L in this Section will denote the Yang-Mills values). In the 
modified action we add a source term A a A a together with the source terms 
J^A atl , as before, in order to derive the Feynman rules. We also have to add a 
gauge fixing term, so I will choose a non-covariant gauge that is well-defined 
and does not depend on A , namely the axial gauge fixing term 

L ax = — — (15) 

with a purely spatial four-vector n M = (0,n). 

The usual Feynman rules for Yang-Mills with action S in the axial gauge 
involve the gauge field propagator 
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and the usual QCD vertices [5]. The new Feynman rules that we get for the 
modified action, using (8), involve an A — A propagator 



Goo — Goo ~ G c 



(17) 



(I will not denote the color indices where obvious) and A — A as well as A — A 
propagators 



Goi and G^ are the same as the usual Yang-Mills. That is, the effect of the 
classical constraint of Gauss's law has been to split the Coulomb interaction 
from the propagator and treat it classically. We also have the usual vertices of 
QCD, and an additional set of vertices: for every QCD vertex that contains 
A we have a vertex where an A leg is replaced by A. There is also the 
ghost sector which is similar to the previous Section, and will be described 
shortly. This ghost sector is not related to the usual Fadeev-Popov ghosts, 
which I will assume that decouple [5]. We can now discuss renormalizability 
and relation to the usual Yang-Mills: 

As far as the gauge field terms A 2 , A 3 and A A terms are concerned we get 
the same terms as usual Yang-Mills: what used to be a Goo propagator can 
now be formed as either the sum of Goo, two Goa and one G\\, or Goo an d 
one Goa, in both cases equal to the usual Goo- This can be verified for every 
diagram. There is one difference: The ghost loop cancels an instantaneous 
Coulomb loop from every diagram in the usual Yang-Mills that contains a 
closed Aq loop. I will assume that regularization can be performed so that 
these loops vanish, that is, relations like 



This assumption is supported by results of split dimensional regularization 
[6]. In fact, even without this assumption, it is possible that the theory will 
be renormalizable, since it is the same closed loop that is missing from every 





(19) 



hold, together with usual axial gauge integrals like 




(20) 
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would-be Yang-Mills diagram, the verification of this, however, would be 
highly non-covariant. In any case, with the previous assumption, the theory 
is renormalizable and equivalent to the usual Yang-Mills as far as the gauge 
fields are concerned. The divergencies of the remaining terms XSS/SAq in 
the modified action are the same as their counterparts in S, as in the previous 
Section, so these are renormalizable too. We now turn to the discussion of 
the A terms. 

Consider first the one-loop terms with external legs A with zero momen- 
tum, for what would be an effective potential term U (A). All the propagators 
can run inside the loop with one exception: The Aq — A and A — A propagators 
cannot appear because the vertices are linear in A. An example of such a di- 
agram is shown in Fig. 1. So the Coulomb interaction is missing. Had these 
terms been there we would have the full covariant propagator G^ v in the 
loop. All these diagrams, therefore, would add up to zero (their value would 
be the same as an effective potential term for A in the usual Yang-Mills 
which does not exist because of gauge invariance). Accordingly the sum of 
the diagrams that contain a Coulomb interaction between two external legs 
A has to be subtracted. The relevant vertices (A — A — Ai) are the same 
as the QCD vertex (A — A — Ai) and the missing terms with a Coulomb 
interaction running between two external legs A correspond to an effective 
interaction term 



(in momentum space, with k the momentum running in the loop) with 



where f acd f bcd = C5 ab , A 2 = A a A a . One can add this effective interaction term 
in the usual Yang-Mills action in order to derive the terms in the effective 
action that depend on A. It corresponds to missing terms so the pieces 
calculated have to be subtracted. The calculation is highly non-covariant, 
however, once we have identified the effect of the missing diagrams as an 
effective interaction term to be added to the usual Yang-Mills action, there 
is no reason to continue in the axial gauge in order to complete the calculation 
[7]. It is more convenient to choose a Feynman gauge-fixing term, (d^A^) 2 /2, 
in which case the effective propagators become: 



k 2 3 



(21) 



m 2 = g 2 C\ 2 



(22) 



-1 



(23) 
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*=H* + w%fe\ (24) 

The usual Fadeev-Popov ghosts do not contribute and the calculation of the 
relevant terms gives the final form of the effective action: 

-Z(X) XD-E- ^Z(X) F 2 + U(X) (25) 
with 

f/(A) = Cl ^A 4 (M^-i) (26) 

Z(X) = l + c 2 g 2 \n^ (27) 
A* 

where c\ = C 2 /647r 2 , C2 = C/6tc 2 . We have an effective action of the 
Coleman- Weinberg form [8] with a few differences: 

First, the potential appears with the opposite sign, since it corresponds to 
missing diagrams. The counterterms that were needed for its renormalization 
were chosen so that U"(0) = and U'(fi) = 0. The condition U'(0) = does 
not fix the counterterms, so we have to pick the scale p, (primes denote 
derivatives with respect to A). 

The factors Z in the first two terms of the effective action (25) can be 
calculated in the same manner, in terms of missing diagrams. They are 
expected to be the same by individual diagram inspection. A first calculation 
gives the result presented above. Z(X) was calculated from the A-dependent, 
p 2 coefficient, of the Aq — A$ wavefunction renormalization diagrams with 
external momentum p (subtracted from the tree level term). Other terms 
(Ai — Aj for example) should give the same value for Z. This, however, 
will have to be verified, because of the asymmetry of the Feynman rules 
described here. In any case, for the preliminary analysis presented below, 
the main fact that we need is that renormalization conditions can be chosen 
so that Z(p) = 1. 

Although there is only one coupling constant in the theory the one-loop 
terms are reliable when g 2 \n(X/p) is small, since the tree value for U is 
and for Z is 1. 

Another possibility that was not encountered in the simple model of the 
previous Section, is the generation of (VA) 2 terms, because of the asymmetry 
of the Feynman rules described here. The tree level value of this term should 
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be set to zero by counterterms, it is a possible fine-tuning problem of the 
method. Generation of these terms at higher order via the Coleman- Weinberg 
mechanism does not affect the analysis presented below. 

At A ~ \i we have Z fa 1 and ordinary perturbative Yang-Mills. As a first 
approximation to the effective action we take the abelian, electric parts, for 
Z «1 

T = / AV 2 A + l-A V 2 A + 17(A). (28) 



The equations 



5£o 
5X 

SA 



= (29) 
= (30) 



have the simultaneous solution 



with \ B the solution of 



A = -Xb (31) 



V*A = f . (32) 

Note, first, that, although similar in form, the two equations (29) and (30) 
are quite different conceptually: The first equation is an expression of the 
classical nature of Gauss's law and should be satisfied identically to all orders. 
The second equation is a usual semiclassical equation, only to be satisfied 
approximately. 

Equation (32) has a soliton (bounce, bubble) solution A#(r), spherically 
symmetric in the three-dimensional radius r, similar to the bounce solutions 
that are associated with tunneling at zero and finite temperature [9, 10], 
although here it is not related to either tunneling or finite temperature. A^(r) 
is of order \i for sufficiently small r, and goes rapidly to zero for r larger than 
the radius Rb of the bounce (which is of order p-)- It corresponds to a 
confining potential for the electric field (31), reminiscent of a bag model, 
such that deep inside the bounce we have ordinary perturbative Yang-Mills 
with a zero electric field, and a strong electric field appearing as we get closer 
to the bounce radius (in fact, the bounce solution is not of the thin wall type, 
it resembles more the three-dimensional thick wall bubbles of [10], and has 
to be determined numerically). 
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As we approach the radius of the bubble, however, when A goes to zero, 
the approximation Z w 1 does not hold; in fact, when A is non-perturbatively 
small Z goes to zero and the full non-perturbative features of QCD be- 
come important [11, 12]. Our perturbative parameter is the same as that for 
Coleman- Weinberg models, namely g 2 ln(A//z), and only when this is small 
can the higher loop effects be neglected. It is, of course, possible to improve 
on these results with the use of the renormalization group, and even when 
this is not the case, the effective action presented here may be quite useful 
phenomenologically. 

Another important fact that should be mentioned for the effective ac- 
tion (25) proposed here, as well as for the "dynamically induced" mass term 
(22), is that they are gauge invariant, provided the auxiliary field A is gauge 
covariant (under a gauge transformation V we have A — > VXV^ 1 ). The conse- 
quences of this gauge invariance, and its possible associated BRS symmetries, 
are not obvious because of the peculiarities of the Feynman rules described 
here; this is, however, another indication in support of the arguments pre- 
sented above, namely that the theory is renormalizable and compatible with 
ordinary, perturbative Yang-Mills. 

4 Comments 

The upshot of this work, as far as the non-abelian gauge theory is concerned 
is the following: there is a sector of the theory, namely the Coulomb interac- 
tion, that is purely classical in nature, by virtue of Gauss's law; this cannot 
be expressed self-consistently in perturbation theory unless one employs a 
skewed set of Feynman rules of the type described here. This has the effect 
of generating additional terms in the effective action that reveal both the ap- 
pearance of the confining properties of the theory and the limits of validity 
of perturbation theory. 

The appearance of the inverted effective potential term U(X) that was 
described in the previous Section does not indicate an energetic instability 
of the theory, it is, however, related to the instability of the perturbative 
vacuum. In fact, there is no direct Hamiltonian interpretation for the effective 
action presented here, A remains an auxiliary Lagrange multiplier that has 
to be eliminated via 5T/5X = 0. The perturbative vacuum X — /i, however, 
cannot exist for all space, since it has infinite action. A finite action soliton 
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solution that was presented in the previous Section shows signs of confining 
behaviour and, at the same time, gives clues about the region of validity of 
perturbation theory. 

It is not clear whether this method hints on fundamental results (problems 
rather) of the usual quantization procedure of non-abelian gauge theories, or 
is purely of a heuristic value, I would like, however, to make some comments 
in support of the former: First of all, something that is obvious, but should 
be, nevertheless, mentioned, is that this method does not change at all the 
rules for the abelian gauge theory. The Coulomb interaction splits again as 
presented here. Since, however, photons do not have any self-interactions, 
the result is equivalent to the usual Feynman rules for abelian gauge theory 
(the inclusion of fermions is also straightforward). The usual Fadeev- Popov 
procedure (abelian or non-abelian) also is not related and does not change 
by the method presented here. However, in the process of quantization of the 
non-abelian gauge theory, while changing from the Hamiltonian formalism 
with non-covariant gauge fixing to the Lagrangian formalism with covariant 
gauge fixing there are many manipulations of constraints such as Gauss's 
law. The treatment of these constraints is not completely justified from the 
present point of view. In many cases the actual, physical, fields are used 
as Lagrange multipliers, and a relation such as (8) does not appear. It is 
possible that some important piece of information of the theory is lost in the 
process. 

5 Discussion 

In this work I developed a formalism in order to treat interactions between 
quantum and classical fields through the path integral. It seems that one is 
able to describe quantum-classical interactions in field theory self-consistently 
with this method. The path integral formalism is particularly simple and can 
hopefully be generalized in other cases involving such interactions. 

Some further applications of this work, besides the conceptual problem 
of quantum-classical interactions, would be in cases of a non-renormalizable 
classical interaction, such as gravity. Note that, in the example described in 
Sec. 2, there is nothing that demands the classical field, 0, to have a renor- 
malizable quantum equivalent. Even a non-renormalizable self-interaction of 
the classical field does not generate higher order terms since there are no 
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classical loops. There may be of course quantum-classical interaction terms 
that generate higher order terms, so the problem of renormalizability then 
depends on the specifics of the model. Even so, this method may be useful 
in treating the stress-tensor renormalization equations 

= -8nG < T M „ > (33) 

in a self-consistent way. 

I have also given an example of an application of this work in the case 
of the non-abelian gauge field theory by using this method in order to treat 
Gauss's law as a classical equation. This is different in spirit than the previous 
discussion, it shows, however, another application of this method in this case, 
where it can be useful for studying the infrared properties of the theory. 
The treatment of the Coulomb interaction as purely classical in nature has 
resulted in the generation of purely quantum terms of the Coleman- Weinberg 
type that are related to the confinement mechanism. 

Even though the treatment of the non-abelian gauge theory presented 
here has a different motivation than previous works, there are several com- 
mon features with various other approaches to confinement [11, 12, 13, 14, 
15, 16, 17, 18]. Namely we find: the appearance of a negative metric propa- 
gator (the A — A line which is crucial in describing the classical nature); the 
emergence of a "dynamical" mass term (22) and an effective potential gen- 
erated through radiative corrections; an effective action reminiscent of older 
phenomenological actions that were used to describe confinement (although 
with important differences described before). 

This method can probably be useful in various other problems in theories 
that have infrared singularities: one can, presumably, use this method in 
order to study the infrared modes of these theories classically, and include 
the quantum corrections from higher momentum modes self-consistently. The 
problem here, of course, is the use of an arbitrary cut-off scale, and care 
should be taken in order to derive cut-off independent results. 



Acknowledgements 

This work was done in the National Technical University of Athens. I am 
grateful to the people of the Physics Department for their hospitality. 



13 



References 

[1] E. Gozzi, M. Reuter and W. D. Thacker, Phys. Rev. D40, 3363 (1989). 

[2] R. Penko and D. Mauro, Eur. J. Phys. 27, 1241 (2006). 

[3] A. H. Mueller and D. T. Son, Phys. Lett. B582, 279 (2004). 

[4] Sangyong Jeon, Phys. Rev. C72, 014907 (2005). 

[5] G. Leibbrandt, Rev. Mod. Phys. 59, 1067 (1987). 

[6] G. Leibbrandt and J. Williams, Nucl. Phys. B475, 469 (1996). 

[7] R. Jackiw, Phys. Rev. D9, 1686 (1974). 

[8] S. Coleman and E. J. Weinberg, Phys. Rev. D7, 1888 (1973). 

[9] S. Coleman, Aspects of Symmetry, Cambridge Univ. Press (1985). 
[10] A. Linde, Nucl. Phys. B216, 421 (1983). 
[11] J. Kogut and L. Susskind, Phys. Rev D9, 3501 (1974). 
[12] M. Bander and P. Thomas, Phys. Rev D12, 1798 (1975). 
[13] G. 't Hooft, Nucl. Phys. Proc. Suppl. 121, 333 (2003). 
[14] D. Zwanziger, Nucl. Phys. B518, 237 (1998). 
[15] J. M. Cornwall, Phys. Rev. D26, 1453 (1982). 
[16] M. Lavelle and D. McMullan, Phys. Rept. 279, 1 (1997). 
[17] P. Watson and R. Alkofer, Phys. Rev. Lett. 86, 5239 (2001). 
[18] C. S. Fischer, J. Phys. G 32, R253 (2006). 



14 



Fig. 1: An example of a diagram used for the calculation of U(X) and the 
derivation of Eqs. (21, 22). The wiggly lines denote A i: Aj fields. The curly 
lines denote the A field and the solid lines denote A . The solid lines with a 
dot denote the modified A — A propagator, Goo; derived in the text, that 
does not contain the Coulomb interaction; the A — A and A — A lines that 
carry the Coulomb propagator cannot be also added in its place, because the 
vertices are only linear in A. 
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